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Abstract 

We study a cosmological model in which phantom dark energy is coupled to dark matter 
by phenomenologically introducing a coupled term to the equations of motion of dark energy 
and dark matter. This term is parameterized by a dimensionless coupling function 6, Hubble 
parameter and the energy density of dark matter, and it describes an energy flow between the 
dark energy and dark matter. We discuss two cases: one is the case where the equation-of-state 
LOe of the dark energy is a constant; the other is that the dimensionless coupling function 5 is a 
constant. We investigate the effect of the interaction on the evolution of the universe, the total 
lifetime of the universe, and the ratio of the period when the universe is in the coincidence state 
to its total lifetime. It turns out that the interaction will produce significant deviation from the 
case without the interaction. 
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1 Introduction 



A lot of evidence of astronomical observations indicates that our universe is currently in ac- 
celerated expansion. This is first revealed by observing high red-shift supernova la Cross 
checks confirm this from the cosmic microwave background radiation and large scale struc- 
ture ini- To explain this accelerated expansion, some proposals have been suggested recently 
in the literature, for instance, by modifying Einstein's general relativity in the cosmic distance 
scale, employing brane world scenario and so on. 

In Einstein's general relativity, however, in order to give an explanation for the accelerated 
expansion, one has to introduce an energy component to the energy density of the universe 
with a large negative pressure, which drives the universe to accelerated expand. This energy 
component is dubbed dark energy in the literature. All astronomical observations indicate that 
our universe is flat and it consists of approximately 72% dark energy, 21% dark matter, 4.5% 
baryon matter and 0.5% others like radiation, etc. A simple candidate of the dark energy is a 
tiny positive cosmological constant, which was introduced by Einstein in 1917, two years later 
since he established general relativity. If the dark energy is the cosmological constant, one has 
to answer the question why the cosmological constant is so small, ~ 10~^^^(Mp)^, rather than 
~ (Mp)^, which is expected from local quantum field theory Here Mp ~ W^^Gev is the 
Planck mass scale. Although the small cosmological constant is consistent with all observational 
data so far, recall that a slow roll scalar field can derive the universe to accelerated expand in 
the inflation models, it is therefore imaginable to use a dynamical field to mimic the behavior of 
the dark energy. In particular, it is hoped that one can solve the so-called coincidence problem 
by employing a dynamical field. The model of scalar field(s) acting as the dark energy is called 
quintessence model [Sj. Following A:-inflation model 6 , it is also natural to use a field with 
noncanonical kinetic term to explain the currently accelerated expansion of the universe. Such 
models are referred to as /c-essence models with some interesting features Suppose that the 
dark energy has the equation of state, pe = uiePe-, where pe and pe are pressure and energy 
density, respectively. In order to derive the universe to accelerated expand, one has to have 
LOf. < —1/3. Note that for the cosmological constant, uJe = —1; for the quintessence model, 
— 1 < We < —1/3; and for the fc-essence model, in general one may has Ue > —1 or uje < —1, but 
it is physically implausible to cross = — 1 [Hi- 
lt is well known that if < —1, the dark energy will violate all energy conditions [5]. 
However, such dark energy models [TO] are still consistent with observation data (—1.46 < 
iOf, < —0.78) ir. The dark energy model with tUe < —1 is called phantom dark energy model. 
One remarkable feature of the phantom model is that the universe will end with a "big rip" 
(future singularity). That is, for a phantom dominated universe, its total lifetime is finite (see 
also jl21ll3j ^. Before the death of the universe, the phantom dark energy will rip apart all bound 
structures like the Milky Way, solar system. Earth, and ultimately the molecules, atoms, nuclei, 
and nucleons of which we are composed [13] . 
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Usually it is assumed that the dark energy is coupled to other matter fields only through 
gravity. Since the first principle is still not available to discuss the nature of dark energy and 
dark matter, it is therefore conceivable to consider possible interaction between the dark energy 
and dark matter. Indeed there exist a lot of literature on this subject (see for example |151 1161 
I17[ I18 | [T!I| and references therein). In this paper, we also consider an interaction model between 
the dark energy and dark matter by phenomenologically introducing an interaction term to the 
equations of motion of dark energy and dark matter, which describes an energy flow between 
the dark energy and dark matter. We restrict ourselves to the case that the dark energy is a 
phantom one. Constraint from supernova type la data on such a coupled dark energy model has 
been investigated very recently JH] (see also |17| I18j). Here we are interested in how such an 
interaction between the phantom dark energy and dark matter affects the evolution and total 
lifetime of the universe. 

On the other hand, one important aspect of dark energy problem is the so-called coincidence 
problem. Roughly specking, the question is why the energy densities of dark energy and dark 
matter are in the same order just now. In other words, we live in a very special epoch when 
the dark energy and dark matter densities are comparable. Most recently, developing the idea 
proposed by Mclnnes [211] . Scherrer [25 has attacked this coincidence problem for a phantom 
dominated universe. Since the total lifetime of the phantom universe is finite, it is therefore 
possible to calculate the fraction of its total lifetime of the universe for which the dark energy 
and dark matter densities are roughly comparable. It has been found that the coincidence 
problem can be significantly ameliorated in such a phantom dominated universe in the sense 
that the fraction of the total lifetime is not negligibly small. In this paper we will also study 
the effect of the interaction on the ratio of the period to its total lifetime when the universe is 
in the coincidence state. 

The organization of this paper is as follows. In the next section, we first introduce the 
coupled dark energy model. In Sec. 3 we discuss the case where the phantom dark energy has 
a constant equation-of-state uJe- In Sec. 4 we study the case with a constant coupling function 
6 introduced in Sec. 2. In this case, the equation-of-state cUe will no longer be a constant. The 
conclusion and discussion are presented in Sec. 5. 

2 Interacting phantom dark energy with dark matter 

Let us consider a csomological model which only contains dark matter and dark energy (gen- 
eralizing to include the baryon matter and radiation is straightforward). A phenomenological 
model of interaction between the dark matter and dark energy is assumed through an energy 
exchange between them. Then the equations of motion of dark matter and dark energy in a flat 
FRW metric with a scale factor a can be written as 
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Pe + 3H{pe + Pe) = -6Hp^, (2.2) 

where pm and Pm are the energy density and pressure of dark matter, while pe and Pe for dark 
energy, H = a/a\s the Hubble parameter, and (5 is a dimensionless coupling function. Suppose 
that the dark matter has pm = and the dark energy has the equation of state pe = oJePe- Note 
that in general is a function of time, rather than a constant. Clearly the total energy density 
of the universe, pt = Pm + Pe, obeys the usual continuity equation 

Pt + 3H{pt+pt) = 0, (2.3) 

with the total pressure pt = Pe- The Friedmann equation is 

= -^Pu (2.4) 
and the acceleration of scale factor is determined by the equation 



{pt + 3pi) , (2.5) 



a 3 

where G is the Newton gravitational constant. 

In general the coupling function 5 may depend on all degrees of freedom of dark matter and 
dark energy. However, if 5 is dependent of the scale factor only, one then can integrate (|2.H) 
and obtain 

Pm = P™,oa-3e/'^'^", (2.6) 

where a = log a and pm,o is an integration constant. Substituting this into (|2.2|) . in order to 
get the relation between the energy density of the dark energy and scale factor, one has to first 
be given the relation of the pressure to energy density of the dark energy, namely the equation- 
of-state (jJe- Here we follow another approach to study the cosmological model by assuming a 
relation between the energy density of dark energy and that of dark matter as follows JH] '■ 

r^P^=P^(^\\ (2.7) 

Pm Pm,0 \«0/ 

where pefl, ao and ^ are three constants. Set the current value of the scale factor be one, namely 
ao = 1, then />e,o and pm,o can be interpreted as the current dark energy density and dark matter 
energy density, respectively. 

In this paper we will consider two special cases. One is the case where uJe is kept as a 
constant. The other is the case where the coupling function 5 is a constant. 

3 Cosmology with a constant equation-of-state of phantom dark 
energy 

In this section we consider the case with a constant LOe- From the relation p.7|) we have 

Pt, Pm = 1 , . . Pt, (3.1) 
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where the constant A = Pe,o/pm,o = ^e,o/^m,o, ^e,o and r^m.o ai"e the fractions of the energy 
densities of dark energy and dark matter at present, respectively. The total energy density 
satisfies 

dpt 3 1 + (1 + uJe)Aa^ 



Integrating this yields 



Pt = Pt,oa''[l - f^e,o(l - a«)]-'"^/^ (3.3) 

where the constant pts) = Pe,o + Pm,o- Therefore the Friedmann equation can be written down 
as 

= Hla-\\ - f)e,o(l - a«)]-3-^/«, (3.4) 

with i?o being the present Hubble parameter. By using of p.lj) and (|3.3() . one can get the 
coupling function b from (|2.H) . 

. g , (^ + 3a;e)^a^ ^ + ZuJe ^^ 
5 = 3 + -J — = T— -TT — = Pe, 3.5) 

where an over dot denotes the derivative with respect to the cosmic time i. This can be expressed 
further as 

SZe,o + (1 - "e,o)a * 

where = — f^e,o(^ + Swg). We have 5(a ^ 1) = 5o and S{a — > oo) = (^o/^e,o- Therefore we 
see that when ^ > — Swg, < 0, which implies that the energy flow is from the dark matter to 
dark energy. On the contrary, when < ^ < — Swe, the energy flow is from the phantom dark 
energy to dark matter. This can also be understood from the equations of motion (|2.H) and 
1)2. 2|) . Further, we can see from 1)3. 5() that there is no coupling between the dark energy and dark 
matter as ^ = — Swg. Of course this is true only for case where We is a constant. In addition, we 
can see from (|3.1|) and (|3.3|) that in this model, the universe is dominated by the dark matter 
at early times, while dominated by the phantom dark energy at later times. 
The deceleration parameter q is 

aa ii 3 1 - rie,o + (1 + '^e)^^e,oa^ 

H'^ 2 1 — i2e^o(l — a^) 

Note that qia — > oo) = —1 + 3(1 +a;e)/(2f^e,o) and q{a ^ 1) = — 1 + 3a;ef^e,o/2, they are always 
negative because < and q[a ^ 1) < q[a oo). In Fig. 1-3 we plot the relation of the 
deceleration parameter to the red shift defined by z = 1/a — 1 for the different lOe and ^. In 
plots we take the fraction of the dark energy r2e,o = 0.72. From figures we can see that for the 
case with a fixed We, a larger ^ leads to a smaller red shift when the universe transits from the 
deceleration phase to acceleration phase. On the other hand, for the case with a fixed ^, a larger 
We has a smaller red shift for that transition from the deceleration to acceleration phase. 
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Figure 1: The deceleration parameter q versus the red shift z for the case of Qe,o = 0.72 and 
We = —1.5. Three curves from top to bottom correspond to the cases ^ = 5.5, 4.5 and 3, 
respectively. Note that the case of ^ = 4.5 = —3uJe is just the case without interaction. 

q 




Figure 2: The deceleration parameter q versus the red shift z for the case of J7e,o = 0.72 and 
We = —1.3. Three curves from top to bottom correspond to the cases ^ = 5.5, 4.5 and 3, 
respectively. 



The total lifetime of the universe can be obtained by integrating the Friedmann equation 
(jS31). It is 

POO 

tu = H,^ / da aV2[i _ o(i _ a«)]3-^/2«. (3.8) 
Jo 

Here we are interested in the change of the lifetime due to the interaction between the dark 
energy and dark matter. Note that when = —SuJe, the interaction disappears. Denote the 
total lifetime by tx for this case, one has 

roo 

tT = H^^ / da ai/2[i _ o(i _ a~^^^)]~^/\ (3.9) 
Jo 

Denote the ratio of the lifetimes tu to tj- by g: 

g ^ t_u_ CdaaV2[i_f]^_o(i_a?)]3-e/2e 



tr ~ /o°° da aV2[i _ ^(i _ a~3..e)]-i/2 

Hr^l^t'C^ - f]e,0)"(^+^=)/2'^^r-V2^e-l(l + r)-l/2^^' 



(3.10) 
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Figure 3: The deceleration parameter q versus the red shift z for the case of Qe,o = 0.72 and 
We = —1.1. Three curves from top to bottom correspond to the cases ^ = 5.5, 4.5 and 3, 
respectively. 



In Fig. 4 we plot the ratio g for three different equation-of-states uJe 



-1.5, -1.3 and -1.1. 




Figure 4: The ratio g of total lifetimes versus the parameter ^ = x for the case of r2e,o = 0.72. 
Three curves from bottom to top correspond to the cases uJe = —1.5, —1.3 and —1.1, respectively. 

Clearly, for a fixed uJe, the universe with a larger ^ has a longer lifetime, while for a fixed ^, a 
larger coe leads to a longer lifetime. Note that in Fig. 4 the three points where three curves cross 
the ^ axis correspond to the situation (.^ = —3uJe) without interaction between the dark energy 
and dark matter. In Fig. 5 the ratio g is plotted versus the parameters ^ and uJe- Note that for 
a constant equation-of-state uJe, the total lifetime of the universe approximately is |21j : 



tT 



l+UJe 



where tm is the age of the universe when the matter and phantom dark energy densities are 
equal. When uJe = —1.5, —1.3 and —1.1, one has t-r = 3tm, 4.3tm; and Htm, respectively. In 
order to satisfy the observation data, ^ has to be chosen at least so that g > 1/3, 1/4.3 and 
1/11, respectively. We can see from Fig. 4 that the constraint of the age of the universe on the 
model is very weak; the parameter ^ can be as small as 1. On the contrary, if one requires that 
the transition of the universe from the deceleration phase to acceleration phase happens around 
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at red shift z < 1, one can see from Fig. 1-3 that one needs to take ^ > 3. Note that the best fit 
of the ACDM model indicates that this transition happens at z ~ 0.5. However, it is aUowed 
that the transition happens at z G (0.3 ~ 1) in the different models. Clearly it is quite necessary 
to make a detailed numerical analysis and to give constraints on the parameters of the coupled 
dark energy model by using of the data from supernova, cosmic microwave background radiation 
and large scale structure. Here let us just mention that a partial analysis of the model according 
to the data of supernova has been done more recently jT5]. The results shows that the supernova 
data favor a negatively coupled phantom dark energy with uje < —1 and — 10 < 5o < — 1- 




Figure 5: The ratio g of total lifetimes versus the parameters ^ = x and Ug in the case of 
Jle,o = 0.72. 

Next we turn to the coincidence problem. Following j21j . we calculate the ratio of the period 
when the universe is in the coincidence state to the total lifetime of the universe. That is, we 
will calculate the quantity 

/ = ^, (3.11) 

tu 

where tc is defined by 

tc = Ho^ daa^/'^[l-ne,oil-a^)f^^^'^^. (3.12) 

J ai 

During the coincidence state, the energy density of dark energy is comparable to that of dark 
matter and the scale factor evolves from ai to 02. What is the exact meaning by the term 
"comparable" ? This is not a well-defined question in order to determine the scale factors ai and 
a2. In [2^, Scherrer defined a scale of the energy density ratio rg so that the dark energy and 
dark matter densities are regarded as to be comparable if they differ by less than the ratio tq in 
either direction. He found that the ratio varies from 1/3 to 1/8 as varies from —1.5 to —1.1 
if ro = 10 in a phantom dark energy model without interaction between the dark energy and 
dark matter. In this sense indeed the coincidence problem is significantly ameliorated in the 
phantom model because the ratio is not so small. As a result it is not so strange that we live in 
the epoch when the energy densities of dark energy and dark matter are of the same order. Now 
we want to see how the fraction varies when the phenomenological interaction is introduced. 
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The fraction of the total hfetime of the universe for which the universe is in the coincidence 
state, turns out to be 



/ 



/i7,,,r3/2«-i(l+r)3-e/2C^^ 



\rl'^^ 2^1 f 3 3.. 1 , 3 



2-fl 



2?' 



3aJe 1 I 
2g ' 2g ' 



-3/25 



2«' 



3i^e 1 I 3_ 

25 ' ^ "I" 25 ' 



--]) 



r[3/2€]r[-3(l+a;e)/25] 

r[-3<^e/25] 



(3.13) 



Note that this ratio is independent of the current density parameter Vte,o- In Fig- 6 and 7 we 
plot the ratio / versus the scale ro for different parameters u)(. and Clearly for the case with 
fixed uJe and ^, a larger ro leads to a larger ratio /. On the other hand, for the case with fixed 
ro and Wg, a smaller ^ gives us a larger ratio. For example, we see from Fig. 6 that when ro = 10 
and uje = —1.5, the ratio / ~ 0.45 for ^ = 3, more large than the case (/ = 1/3) without the 
interaction. Note that the middle curve in Fig. 6 corresponds to the case without the interaction 
= —2>uJe)- In Fig. 8, we plot the ratio / versus the parameters ^ and ro for a fixed uJe = —1.3. 
From Fig. 6-8, one can see that indeed the period when the universe is in the coincidence state is 
comparable to its total lifetime. In addition. We note from (|2.7|) that a larger ^ means that the 
universe will be dominated more quickly by the phantom dark energy, it corresponds to have a 
longer total lifetime of the universe, which can be seen from Fig. 4 and 5. Thus it is natural to 
have a smaller fraction / than the case with a smaller ^. 




Figure 6: The ratio / versus the parameter ro for the case of oJe = —1.5. Three curves from top 
to bottom correspond to the cases of ^ = 3, 4.5 and 5.5, respectively. 



4 Cosmology with a constant coupling parameter 

In this section we consider the case with a constant coupling function 6. In this case, we have 
the energy density of dark matter 

Pm = Pm,oa~^^^. (4.1) 
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Figure 7: The ratio / versus the parameter ro for the case of oJe = —1.1. Three curves from top 
to bottom correspond to the cases of ^ = 2, 3.3 and 4, respectively. 




Figure 8: The ratio / versus the parameters tq and ^ for the case of uje = —1.3. In the figure 

And then the dark energy density has the relation to the scale factor 

Pe = Pefla-''^^^^. (4.2) 

The Friedmann equation turns out to be 

= i?o'(f^m,oa"'+^ + 0^,00-=^+^+^). (4.3) 

In this case, the equation-of-state of the dark energy will depend on time (scale factor). From 
H2.2|) . we can obtain 

When (5 = 0, one has ^ = — 3a'e- This situation is just the case without the interaction. From 
(|4.2|) and (|4.4|) . one can see that in order the dark energy to be phantom, 5 + > 3 has to 
be satisfied so that the dark energy density increases with the scale factor. When 5 + ^ = 3, 
although the dark energy density keeps as a constant, it does not act as a cosmological constant 
due to the interaction between the dark energy and dark matter. We see from ()4.4p that 

LOe{a ^ 0) = -sign{d) ■ 00, uJe{a ^ I) = 5 ^Ti — ' LOe{a ^ 00) = — . (4.5) 

o o s ig "J 
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It is interesting to note that uje diverges as a ^ 0. This is due to the fact that the dark energy 
density (|4.2)) approaches to zero very quickly, while the energy density of dark matter (|4.H) goes 
to infinity as a ^ 0, in order for the equation H2.2|) to hold, the parameter has to go to 
infinity. In fact, the dark energy does not play any role at early times, which can be seen from 
the Friedmann equation (|4.3j) . Therefore here the divergence of does not make any sense in 
physics. 

The deceleration parameter is found to be 

1 (3 - 5)Vtmfl + (3 - <5 - mefla^ 



1 = -1 + 

5 — ^ile,o)/2 when a = 1 and q - 



(4.6) 



which has q = —1 + {3 — 6 — ^ile,o)/2 when a = 1 and q = {1 — 6 — (,)/2 when a oo. In 
Fig. 9 and 10 we plot the deceleration parameter versus the red shift for the case with a fixed 
^ = 4, different 5, and the case with a fixed 5 = 0.3, different ^, respectively. Note that the 
coupling parameter 6 will affect the expansion law (|4.1|) of dark matter density. We do not 
expect that the usual behavior {pm ~ a^^) will be changed much. So we take the value of 6 so 
that the exponent is changed within 10%. That is, 6 is taken to be in the range of (—0.3,0.3). 
Of course, in order to get a correct constraint on the parameter 6 from the observation data, a 
detailed numerical analysis has to be done. We see from Fig. 9 that for the case with a given 
^, a larger 5 leads to a smaller red shift when the universe transits from a deceleration phase to 
an acceleration phase, while Fig. 10 tells us that for the case with a fixed 5, a larger ^ gives us 
a smaller red shift. 




Figure 9: The deceleration parameter q versus the red shift z for the case of il.e,o = 0.72 and 
= 4:. Three curves from top to bottom correspond to the cases 6 = 0.3, and —0.3, respectively. 



From (|4.3)1 we can get the total lifetime of the universe 



tu = Hq 



-1 







da a ^(r2m,oa '^^^ + ^efiO- 



-3+<5+€<j-l/2 



(4.7) 



We now consider the effect of the interaction on the total lifetime. Note that the total lifetime 
of the universe without the interaction is 

coo 



poo 

tT = Hq^ / da a-^{Qrn oa~^ + oa-^+^)-^/^. 
Jo 



(4.8) 
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Figure 10: The deceleration parameter q versus the red shift z for the case of Oe.o = 0.72 and 

5 = 0.3. Three curves from bottom to top at the q axis correspond to the cases ^ = 5, 4 and 3, 
respectively. 

Denote the ratio tjj/tT by g, we can express this as 

/- dr(i^)(3-'5)/2C^(3-<5-20/2C(i + ^)-l/2 
d^(i:^)3/2C^(3-20/2€(l + ^)-l/2 

In Fig. 11, the ratio g is plotted versus the parameters ^ and 5 for the case Qe,o = 0.72. We see 
that the case (5 > is quite different from the case of 5 < 0. For the case with a fixed 5 > 0, 
a larger ^ leads to a longer lifetime of the universe. On the contrary, for the case with a fixed 

6 < 0, a smaller ^ gives us a longer lifetime. Further, for the case with a fixed ^, a smaller 6 has 
a longer lifetime of the universe. 



Figure 11: The ratio g of total lifetimes versus the parameters ^ = x and (5 = d for the case of 
Oe,o = 0.72. 

Finally we consider the fraction of the total lifetime for which the universe is in the coinci- 
dence state. As the case with a constant uJe considered in the previous section, we calculate the 
following ratio: 

ro ^^^{3-5-20/2?^ +^)-l/2 

'' /~dr r(3-'5-20/2«(i + r)"i/2 • ^^■^''^ 



12 



Figure 12: The ratio / versus the parameter tq for a fixed = 4. Three curves from top to 
bottom correspond to the cases of (5 = 0.3, and —0.3, and corresponding equation-of-state are 
ujefi = —1-47, —1.33 and —1.19, respectively. 

In Fig. 12 we plot the ratio / versus the scale ro for a fixed ^, but different 6. It shows that 
a larger 5 gives a larger ratio for a fixed tq. On the other hand, we plot the ratio / in Fig. 13 
versus the scale rg for a fixed 5, but different ^, which shows that a larger ^ gives us a larger 
ratio for a fixed rg. Note that the case with tUefl = -1-80 is already ruled out |llj . Here we draw 
this case just for illustration. Fig. 14 shows the relation of the ratio / versus the parameters ^ 
and 6 for a fixed scale vq = 5. 

f 

0.i[ 



■ 6 8 10 " 

Figure 13: The ratio / versus the parameter rg for a fixed 6 = 0.3. Three curves from top 
to bottom correspond to the cases of ^ = 5, 4 and 3, and corresponding equation-of-state are 
uje,o = -1-80, —1.47 and —1.13, respectively. 



5 Conclusion 

In summary we discuss a cosmological model in which phantom dark energy has an interaction 
with dark matter. The interaction is introduced phenomenologically by considered an additional 
term [see (|2.1|) and (|2.2j) ] in the equations of motion of dark energy and dark matter. This term 
is parameterized by a product of a dimensionless coupling function 5, Hubble parameter and the 
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Figure 14: The ratio / versus the parameters 5 = d and ^ = x for a given ro = 5. 



energy density of dark matter, and it describes an energy flow between the dark energy and dark 
matter. We discuss two cases, one is the case where the equation-of-state uje = Pe/Pe of the dark 
energy is kept as a constant; the other corresponds to the case with a constant couphng function 
6. We investigate the effect of the interaction on the evolution of the universe, the total lifetime 
of the universe, and the fraction of the total lifetime of the universe for which the universe is 
in the coincidence state, where the energy densities of the dark energy and dark matter are 
comparable. We find that the interaction has rich and significant consequences on these issues. 
For example, the fraction of the total lifetime of the universe for which the universe is in the 
coincidence state can approximately reach 0.45 if we take iVg = —1.5, tq = 10 and ^ = 3. 
It means that the period when the energy densities of the dark energy and dark matter are 
comparable is significantly long, compared to its total lifetime. Thus, it is not so strange that 
we now live in the coincidence state of the universe. In this sense the coincidence problem can 
indeed be significantly ameliorated in such an interacting phantom dark energy model. Finally 
let us stress that the constraints on the parameters of the coupled dark energy model from the 
supernova la data has been analyzed recently ^Hl (see also jElE]), the values of parameters 
we take in this paper are all in the allowed region. Certainly it is not enough to consider the 
constraints from the supernova data only. We expect that the data from the cosmic microwave 
background radiation and in particular, from the large scale structure will give more restrictive 
constraints on this coupled dark energy model. This issue is currently under investigation. 

Note added. After the manuscript was put on the net, we have been informed that a similar 
idea has also been considered by Z.K. Guo and Y.Z. Zhang. In addition, we have noticed that 
following 12^, the coincidence problem has been discussed in a scalar field dark energy model 
with a linear effective potential |22j . 
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